Realizing Colloidal Artificial Ice on Arrays of Optical Traps 
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We demonstrate how a colloidal version of artificial ice can be realized on optical trap lattices. 
Using numerical simulations, we show that this system obeys the ice rules and that for strong colloid- 
colloid interactions, an ordered ground state appears. We show that the ice rule ordering can occur 
for systems with as few as twenty-four traps and that the ordering transition can be observed at 
constant temperature by varying the barrier strength of the traps. 
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In certain spin models, the geometric spin arrange- 
ments frustrate the system since not all of the nearest 
neighbor spin interaction energies can be minimized si- 
multaneously JjL A classic example of this is the spin 
ice system |2, y|, named after the similarity between 
magnetic ordering on a pyrochlore lattice and proton or- 
dering in water ice Spin ice behavior has been ob- 
served in magnetic materials such as Ho2Ti207, where 
the magnetic rare-earth ions form a lattice of corner- 
sharing tetrahedra 0. The spin-spin interaction energy 
in such a system can be minimized locally when two spins 
in each tedrahedron point inward and two point outward, 
leading to exotic disordered states 0. There are several 
open issues in these systems, such as whether long range 
interactions order the system, or whether the true ground 
state of spin ice is ordered Q . 

In atomic spin systems, the size scale is too small to ex- 
amine ordering on the individual spin level directly, and 
very low temperatures are required to freeze the spins. 
Artificial versions of spin ice systems that overcome these 
limitations would be very useful. In a recent experiment, 
a geometrically frustrated system was constructed from 
a square lattice of small, single-domain magnetic islands 
. Each vertex of the lattice represents a meeting point 
for four spins. Wang et al. demonstrated that the sys- 
tem obeys the "ice rules" of two-spins-in, two-spins-out 
at each vertex for closely spaced islands, and has a ran- 
dom spin arrangement for widely spaced islands. Unlike 
in atomic systems, it is possible to image the ground state 
of the resulting spin ice directly using a scanning probe. 

Here, we propose another version of an artificial spin 
ice system in which both statics and dynamics can be 
probed directly. We use numerical simulations to show 
that square ice as well as other frustrated states can be 
constructed using interacting colloidal particles confined 
in two-dimensional (2D) periodic optical trap arrays. 
Due to the micron size scale of the colloids, the ordering 
and dynamics could be imaged with video-microscopy in 
an experiment. The colloidal system may also equilibrate 
much more rapidly than the nanomagnet system, since 
thermal fluctuations are present and can be controlled 
by changing the relative strength of the optical traps. In 
addition, the colloidal interaction can be changed from 



nearest neighbor to longer range simply by adjusting the 
screening length. A variety of different static and dy- 
namical trap geometries can be constructed with optical 
arrays jSj] , and colloidal crystallization and melting have 
already been demonstrated in square and triangular op- 
tical trap arrays |3, 0| . It is also possible to make ar- 
rays with elongated traps that have a double well shape 
such that a single colloid can be located in either well 
[llllT3 |. Colloid-colloid interaction forces between neigh- 
boring traps were strong enough to induce a zig-zag or- 
dering and permit signal propagation in an exp eriment 
on a chain of 23 double well elongated traps jTj . 

We consider an artificial ice system created with 
charged colloids on a square lattice of elongated optical 
traps at a one-to-one filling. Each trap has a double well 
potential similar to those created experimentally 0, ^| 
The vertices where four traps meet correspond to the 
oxygen atoms or the pyrochlore tetrahedrons, and the 
charged colloidal particles model the 'in' or 'out' spins. 
We classify the resulting six vertex types 14] according 
to their electrostatic energy, and study the change in the 
occupancy of different vertex types as a function of col- 
loid charge and trap spacing. For noninteracting colloids 
we find random occupancy of each vertex type. When the 
colloid charge is increased, ice-rule obeying vertices dom- 
inate the system. At high interaction strengths, we find 
a long-range ordered minimum energy state predicted 
previously for spin ice systems j6[. We can recover the 
random vertex occupancy at high colloid charge by in- 
creasing the spacing between traps. When the traps are 
gradually biased in one direction, we obtain a transition 
between two ice-rule obeying ground states. Since it can 
be difficult to construct very large arrays of traps exper- 
imentally, we show that a system with as few as 24 traps 
and open boundary conditions still exhibits the ordering 
transition observed for larger arrays. Previous experi- 
ments focused on order-disorder transitions for colloids 
as a function of increasing trap strength [lol |. Thus, we 
also demonstrate that the spin ordering transition occurs 
at constant temperature when the trap barrier strength 
is varied. 

We perform 2D Brownian dynamics (BD) simulations 
for systems of two sizes. System A contains N = 1800 
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FIG. 1; (a) Schematic of the basic unit cell with four double 
well traps each capturing one colloid, (b-d) Images of a small 
portion of system A with A'' — 1800. Dark circles: colloids; 
ellipses: traps, (b) Random vertex distribution at g = 0. (c) 
Long-range ordered square ice ground state ai q = 1.3. (d) 
Biased system at g = 0.4 with F'^" = 0.02. 



interacting colloids and N = 1800 optical traps with peri- 
odic boundary conditions in the x and y directions. Sys- 
tem B has N = 2A colloids and N — 2A optical traps with 
open boundary conditions. In each case the overdamped 
equation of motion for colloid i is: 



V- 



dt 



(1) 



where the damping constant rj — 1.0. We de- 
fine the unit of distance in the simulation to be uq. 
The colloid-colloid interaction force has a Yukawa or 
screened Coulomb form, F^'^ = "-^o?^ X^i^j 
with V{rij) = (l/ry)exp(-Kry)fy. Here Uj = jr^ -r^l, 
Tij — (r; — rj)/rij, ri(j) is the position of particle 
Fq = Z*^/(47reeo), Z* is the unit of charge, e is the sol- 
vent dielectric constant, q is the dimensionless colloid 
charge, and 1/k is the screening length, where k = 4/ao 
unless otherwise mentioned. We neglect hydrodynamic 
interactions between colloids, which is a reasonable as- 
sumption for charged particles in the low volume frac- 
tion limit. The thermal force is modeled as ran- 
dom Langevin kicks with the properties (F^) = and 
(F'^(t)F^(t')) = 2-qkBT5{t ~t'). Unless otherwise men- 
tioned, = |F-^| = 0. F^^* represents an externally 
applied drive which is set to zero except for the biased 
system, where F^^* = F^'^{yi + y). 

The substrate force F| arises from elongated traps, 
shown schematically in Fig. 1(a), arranged in square 
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0000 


0.001 


IV 


1001 


7.02 


II 


0001 


0.0214 
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1101 


14.977 


III 


0101 


1.0 


IV 


nil 


29.913 



TABLE I: Electrostatic energy Ei/Em for each vertex type. 
An example configuration for each vertex is listed; 1 (0) indi- 
cates a colloid close to (far from) the vertex. 



structures with lattice constant d, as in Fig. 1(b). Each 
trap is composed of two half-parabolic wells of strength 
fp and radius Tp separated by an elongated region of 
length 21 which confines the colloid perpendicular to the 
trap axis and has a small repulsive potential or barrier 
of strength parallel to the axis which pushes the col- 
loid out of the middle of the trap into one of the ends: 

I,. Herer± = |r, - ± Zp|||, 
(r^) is the position of colloid 



p1|| 



{frll){l-ri)Q{l 

i (trap k), and p|| (p'[) is a unit vector parallel (per- 
pendicular) to the axis of trap k. We take 21 = 2ao, 
Tp — 0.4ao, and d ~ 3ao unless otherwise noted. Elon- 
gated traps of this form have been created in previous 
experimental work 0, . Our dimensionless units can 
be converted to physical units for a particular system. 
For example, when ao = 2 //m, e = 2, and Z* — 300e, 
such as in Ref. ISj, Fq = 2.5 pN and the trap ends are 
0.2 /im apart at c? = 3. We find the ground state of each 
configuration using simulated annealing. 

The vertices are categorized into six types, listed in 
Table I, and we identify the percentage occupancy Ni/N 
and energy Ei of each type. Type III and type IV vertices 
each obey the ice rule of a two-in two-out configuration, 
represented here by two colloids close to the vertex and 
two far from the vertex. Locally, the system would prefer 
type I vertices, but such vertices must be compensated by 
highly unfavorable type VI vertices. The colloidal spin 
ice realization differs from the magnetic system, where 
north-north and south-south magnetic interactions at a 
vertex have equal energy. For the colloids, interactions 
between two filled trap ends raise the vertex energy 
whereas two adjacent empty trap ends decrease Ei. Since 
particle number must be conserved, creating empty trap 
ends at one vertex increases the particle load at neigh- 
boring vertices. As a result, the ice rules still apply to 
our system, but they arise due to collective effects rather 
than from a local energy minimization. 

In Fig. 1(b) we illustrate a small part of system A with 
noninteracting colloids at charge q — Q. The distribution 
of Ni/N is consistent with a random arrangement. When 
we increase g to g = 1.3 so that the colloids are strongly 
interacting, we find a nonrandom configuration where the 
system is filled with type III vertices in a checkerboard 
pattern corresponding to the square ice ground state |3 j 
illustrated in Fig. 1(c). We find similar behavior in sys- 
tem B containing only 24 traps, which may be easier to 
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FIG. 2: O: Ni/N; □ : Nn/N; ♦; Nm/N- A: Niv/N; ■: 
Nv/N; •: Nvi/N. (a) iVi/iV vs g at d = 3 and k = 4.0. (b) 
Ni/N vs d at g = 1.3 and k = 4.0. Inset: schematic spin 
representation of the 6 vertex types, (c) Ni/N vs n a,t d = 3 
and q — 1.0. (d) Ni/N vs q for a biased system at d = 3, 
K = 4.0, and F"^" = 0.02. 



realize in an experiment. 

We use Ni/N to map the transition between the ran- 
dom state and the long-range ordered state as a func- 
tion of colloid charge q at fixed trap spacing d — 3 and 
K = 4.0 in Fig. 2(a), as a function of d for fixed q = 1.3 
and K — 4.0 in Fig. 2(b), and as a function of k for fixed 
q — 1.0 and d = 3 in Fig. 2(c). Changing q, d oi k 
changes the relative colloid-colloid interaction strength. 
In Fig. 2(a), Nj/N and Nyi/N decrease with increas- 
ing q for q > 0.1, since this is the most energetically 
unfavorable vertex combination. As q increases above 
q > 0.3, Nji/N, Nv/N and Niv begin to decrease since 
type II, V and IV vertices are also energetically unfavor- 
able. Even though type IV vertices obey ice rules, they 
disappear at high q since in square ice they have higher 
energy than the type III vertices, which ultimately form 
the ground state. 

In Fig. 2(b) we fix g = 1.3 and k, = 4.0 and show that 
changing the trap spacing d from d = 3 to d = 5 produces 
the same transition as changing the colloid charge. As d 
increases, the colloid-colloid interaction strength drops, 
and the vertices that were suppressed by the high colloid 
charge q reappear in the same order: types II, V and IV 
first, followed by types I and VI. If q and d are fixed and 
the inverse screening length n is varied, we find that a 
similar transition from an ordered to a random configu- 
ration occurs, as illustrated in Fig. 2(c) for d = 3 and 
q = 1.0. In experiments with dynamic optical traps, it 
would be straightforward to adjust the trap spacing d, 
allowing both random and ordered limits to be accessed 
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FIG. 3: (a,b) Images of the entire sample for system B, with 
= 24 and open boundary conditions at d = 3, g = 1, 
K = 4.0, and = 0.9. (a) Disordered system at fr = 0.01. 
(b) Square ice ground state at fr = 0.10. (c) The number 
of defects A^^ versus temperature F'^ for the same system 
with fr = 0.03 (filled circles), fr = 0.07 (open squares), fr = 
0.1 (filled diamonds), and /,■ — 0.15 (open triangles), (d) A 
phase diagram showing the ordered and disordered states as 
a function of F^ vs fr. 

readily. The appearance of an ordered state for small 
trap spacing is similar to experimental observations of 
spin liquid ordering under pressure [l6|. 

From these results we see that a colloidal system can 
serve as a model for artificial ice. The advantage of the 
colloid realization is that it is possible to directly observe 
the dynamical behavior of the individual "spins." It is 
also possible to create more complex systems than can 
be achieved with the magnetic system. For instance, to 
create biased traps, we introduce a small driving force 
oriented at a 45° angle along the diagonal of the square 
plaquettes. This causes the colloids to favor sitting 
in the topmost and rightmost halves of the traps, and 
breaks symmetry in the same way as applying a mag- 
netic field alon g th e [100] direction in a pyrochlore sys- 
tem [T7IITIIT9II20II. In Fig. 1(d) we illustrate the ground 
state of a system biased in this way with q = 0.4, d = 3, 
K = 4, and F'^'^ — 0.02. All of the traps are effectively 
tilted in the rightward and upward directions, and the 
system has a ground state made up entirely of type IV 
vertices. In Fig. 2(d), we see that in the biased system, 
increasing q produces a transition between the tilted and 
the spin ice ground states. At the transition, we find 
vertex clusters composed entirely of type III or type IV 
vertices, with pairs of type II and type V vertices located 
at the boundaries between the clusters. Type I and type 
VI vertices never appear in this system. 

In experimental systems using current equipment, a 
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limited number of optical traps are available. For ex- 
ample, 23 traps were employed in recent experiments on 
double well elongated traps . It is therefore pertinent 
to determine how small of a system can still exhibit the 
artificial ice behavior. Additionally, experimental sys- 
tems have no periodic boundary conditions, so it is im- 
portant to understand whether the boundaries affect the 
response of the system. 

In Fig. 3(a,b) we illustrate two configurations of sys- 
tem B, which contains 24 traps and has open boundary 
conditions, for the case q = 1.0, d = 3.0, and k = 4.0 
at a temperature F"^ — 0.9. A disordered configuration, 
shown in Fig. 3(a), results when the strength of the repul- 
sive barrier at the center of each trap, fr = 0.01, is small. 
Here the thermal effects are strong enough that the col- 
loids are hopping between the two different wells in each 
trap. As fr is increased, a freezing transition occurs and 
the system forms the ordered phase, seen in Fig. 3(b) for 
fr = 0.10. To quantify the ordering transition, we mea- 
sure the time-averaged number of defects in the system. 
Nil, which is determined by comparing the colloid config- 
uration C with the two possible ground states for the 24 
traps: G, shown in Fig. 3(b), and G", obtained by flip- 
ping each colloid in Fig. 3(b) to the opposite side of each 
trap. We take = (min(Ef \a - G,|, - G^|)), 

where \Ci — Cf' | = (1) if colloid i is at the same (oppo- 
site) end of the trap in the two configurations. A plot of 
versus in Fig. 3(c) shows that as is decreased, 
there is a freezing transition into the ordered state where 
TVd = 0. The freezing temperature decreases as fr is low- 
ered. The most straightforward experiment to perform 
would be to fix temperature and the other parameters 
while varying the strength of the barrier at the center of 
the traps, fr- To illustrate that the order-disorder arti- 
ficial ice transition can also occur as function of barrier 



strength, we map out a phase diagram as a function of 
fr and F'^ in Fig. 3(d). This shows that for systems as 
small as 24 traps, the artificial ice behavior should be 
experimentally observable. 

All of the systems we have studied up to this point have 
a well defined, long-range ordered ground state. We have 
also found that if a honeycomb arrangement of traps is 
used instead of the square trap arrangement considered 
here, only disordered ground states occur 2lJ. 

In conclusion, we have shown that an artificial ice 
model system can be created using charged colloidal par- 
ticles in arrays of elongated optical traps. The system 
obeys the ice rules and shows a transition between a ran- 
dom configuration and a long-range ordered ground state 
as a function of colloid charge, trap size, and screening 
length. We demonstrate that a thermally induced order- 
disorder transition also occurs in samples with only 24 
traps and open boundary conditions, which should be 
well within the range of current experimental capabili- 
ties. This transition can be observed at fixed temperature 
by varying the trap barrier strength, which would be the 
most straightforward experiment to conduct. Besides op- 
tical traps, other systems including electrophoretic traps 
I22II or patterned surfaces may also be used to confine the 
colloids. Similar effects should occur for vortices in type- 
II superconductors interacting with elongated arrays of 
blind holes. Experimental versions of frustrated colloidal 
systems could allow for direct visualization of the dy- 
namics associated with frustrated spin systems, such as 
deconfined or confined spin arrangements, as well as spin 
dynamics at melting transitions. 

Note added- A simulation study 23] of the experimen- 
tal dipolar system appeared after submission of this work. 

This work was supported by the US Department of 
Energy under Contract No. W-7405-ENG-36. 
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